arXiv:1507.03837vl [physics.optics] 14Jul2015 


Soliton dynamics in gas-filled hollow-core photonic crystal fibers 


Mohammed F. Saleh^’^ and Fabio Biancalana^ 

^School of Engineering and Physieal Seienees, Heriot-Watt University, EHlf 4^^ Edinburgh, UK 
^Department of Mathematies and Physies Engineering, Alexandria University, Alexandria, Egypt 

Gas-filled hollow-core photonic crystal fibers offer unprecedented opportunities to observe novel 
nonlinear phenomena. The various properties of gases that can be used to fill these fibers give addi¬ 
tional degrees of freedom for investigating nonlinear pulse propagation in a wide range of different 
media. In this review, we will consider some of the the new nonlinear interactions that have been 
discovered in recent years, in particular those which are based on soliton dynamics. 


I. INTRODUCTION 

A soliton is a nonlinear localized wave possessing a 
particle-like nature, which maintains its shape during 
propagation, even after an elastic collision with another 
soliton. In optics, this special wave-packet can arise due 
to the balance between nonlinear and dispersive effects. 
Based on confinement in time or space domain, one can 
have either temporal or spatial solitons. Both kinds of 
optical solitons can occur due to the third-order nonlin¬ 
earity, named as Kerr effect [T] that leads to an intensity- 
dependent refractive index of the medium. This nonlin¬ 
ear dependence results in spatial self-focusing and tempo¬ 
ral self-phase modulation [2]. A spatial soliton is formed 
when the self-focusing effects counteracts the natural 
diffraction-induced broadening of the pulse. Similarly, a 
temporal soliton is developed when the self-phase mod¬ 
ulation effects compensates the usual dispersion-induced 
broadening. 

The possibility of soliton propagation in the anomalous 
dispersion regime of an optical fiber was predicted by an¬ 
alyzing theoretically the nonlinear Schrodinger equation 
(NLSE) [3]. Due to the lack of ultrashort pulses in this 
regime, the experimental verification was delayed until 
1980, when Mollenauer et al. were able not only to ex¬ 
cite a fundamental soliton but also a higher-order soliton 

[4] . A Ath-order soliton that represents a joint state of 
N fundamental solitons, can propagate in a periodic way 
consisting of pulse splitting followed by a recovery to the 
original pulse after a characteristic propagation length. 

Supercontinuum generation, which is a massive pulse 
broadening, was first observed in step-index silica-core 
optical fibers by pumping in the normal dispersion regime 

[5] due to mutual interaction between the self-phase mod¬ 
ulation effect and Raman-scattering process [6]. After 
the availability of ultrashort sources in the anomalous 
dispersion regime, experiments show that soliton emis¬ 
sion leads to broadband supercontinuum generation in 
optical fibers I3CI]. In this case, an energetic pulse or a 
Ath-order soliton is continuously temporally compressed 
and spectrally broadened due to the interplay between 
the anomalous dispersion and Kerr effect. After a certain 
propagation distance, the pulse breaks up and a series 
of fundamental solitons are ejected in a process known 
as soliton fission due to higher-order dispersion effects. 
These latter effects are also responsible to stimulate the 


transfer of energy from each fundamental soliton to a 
weak narrowband dispersive wave in the normal disper¬ 
sion regime [12]. The energies of the higher-order soli¬ 
ton was shown to be equal to the sum of energies of its 
constituents m- Because of the intrapulse Raman scat¬ 
tering of silica, the central frequency of the fundamental 
solitons are continuously downshifted during propagation 
IHlIS]. This frequency shift continues until it saturates 
when the soliton gets chirped m- 

Photonic crystal fibers (PCFs) are microstructured 
fibers that can be engineered in various ways to tailor 
their linear and nonlinear properties. PCF structures are 
fabricated using different techniques. The most widely 
used method is based on stacking a number of capillary 
tubes in a suitably shaped preform to form the desired 
arrangement, which is then drawn down to a fiber. This 
has led to the fabrication of the first silica-air endlessly 
single mode PCF m- Followed by a series of successful 
experiments, PCFs with special characteristics have been 
demonstrated such as PCFs with large mode area m, 
dispersion controlled [Hun], hollow-core [21], multicore 
[22], and birefringence [23]. Another different technique 
to manufacture PCFs (so-called Omniguide fibers) is to 
deposit chalcogenide glass on a polymer using thermal 
evaporation, then wrap it to form multilayered Bragg 
fibers [24ll26] . A full thorough review of the history, fab¬ 
rication, theory, numerical modeling, and applications of 
PCFs are found in m- 

Solid-core PCFs guides light via total internal reflec¬ 
tion similar to step-index fibers. However, the additional 
degrees of freedom provided by modifying the capillary 
diameters open up different possibilities to engineer the 
fiber properties, such as shifting its zero dispersion wave¬ 
length (ZDW) [19], or enhancing its Kerr nonlinearity via 
reducing its effective-core area [28]. A number of simul¬ 
taneous experiments have succeeded in exploiting these 
new advantages offered by solid-core PCFs in generating 
broadband supercontinuum [2Ql I291I3T] . Supercontinuum 
generation in solid-core PCFs has been reviewed in [32] 

Hollow-core (HC) PCFs have also attracted much in¬ 
terest since their invention [21] |33j [34] , because of their 
potential for lossless and distortion-free transmission, 
particle trapping, optical sensing, and novel applications 
in nonlinear optics [35H37] . HC-PCFs can be classified 
into two different categories based on the guiding mech¬ 
anism. Photonic bandgap (PBG) HC-PCFs guide light 
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using the concept of forbidden bands in periodic solid 
crystals. If the optical frequency lies within the photonic 
bandgap, propagation through the fiber cladding is pro¬ 
hibited and light is confined inside the hollow-core. These 
fibers offer very low-loss over a restricted wavelength 
bands [38]. As the number of cladding layers increases, 
light confinement becomes stronger. The other category 
is the anti-resonant fibers, such as the Kagome-lattice 
fiber, which can allow waveguiding in a thick low-index 
core surrounded by a thin high-index cladding via the 
anti-resonant Fabry-Perot mechanism [39]. The fibers 
have a broad transmission band in comparison to PEG 
fibers, however, with much higher but tolerable losses. 
Adding more cladding layers does not play a significant 
role in reducing these losses. A lot of work has been re¬ 
cently dedicated towards a dramatic suppression of the 
attenuation coefficient of the anti-resonant fibers, while 
maintaining its broad transmission bandwidth, by induc¬ 
ing fiber-cores with negative-curvatures [MIST]. 

Gas-filled HC-PCFs with Kagome lattice have become 
a powerful alternative for solid-core PCFs especially for 
nonlinear applications [42]. The low-loss wide transmis¬ 
sion window, the pressure-tunability dispersion, and the 
variety of different gases with special properties have of¬ 
fered several opportunities for demonstrating new non¬ 
linear applications [43] such as Stokes generation with a 
drastical reduction in Raman threshold Hi. high har¬ 
monic generation [45] , efficient deep-ultraviolet radiation 
[46] . ionization-induced soliton self-frequency blueshift 
[471150] . strong asymmetrical self-phase modulation, uni¬ 
versal modulational instability m, parity-time symme¬ 
try [52], and temporal crystals ISlISll- The subject of 
this review is to present in detail some of these applica¬ 
tions, especially those which involve soliton dynamics. 

This review paper is organized as follows. In Sec. II, 
we give an overview of the governing equation of non¬ 
linear pulse propagation in Kerr media and the effect 
of Raman-nonlinearity in silica-glass. Secs. Ill and IV 
are dedicated to the study of soliton dynamics in HC- 
PCFs filled by Raman-active and Raman-inactive gases, 
respectively. Conclusions and final remarks are presented 
in Sec. V. 


II. NONLINEAR PULSE PROPAGATION IN 
GUIDED KERR MEDIA 

Nonlinear pulse propagation in a lossless Kerr medium 
can be described by the scalar nonlinear Schrodinger 
equation (NLSE), 

E jm 

—/3„ar^ + 7|^|'A = 0, (1) 

ml 

m=2 

where the slowly varying envelope approximation 
(SVEA) is assumed, A (z, t) is the complex envelope of 
the electric field in units of z is the longitudi¬ 

nal coordinate along the fiber, t is time in a reference 


frame moving with the pulse group velocity, firn is the 
mth order dispersion coefficient calculated at the pulse 
central frequency cjo, and 7 is the nonlinear Kerr coeffi¬ 
cient. This equation is usually integrated using the split- 
step Eourier method m- The second term is usually 
taken as a fit of the dispersion in the frequency domain 00 
as X]m =2 {oj — /m\ or in an approximation-free 

manner as (3 (uj) — fio — fii {uj — ujq) in those cases when 
the dispersion relation /3{uj) of the waveguide under con¬ 
sideration is known. In a regime of deep anomalous dis¬ 
persion, i.e. for 1 /^ 21 ^ \fim> 2 \^ the normalized solution 
of Eq. 0 is the fundamental Schrodinger soliton, 

t/, (^, r) = A^sech {Nt) exp {iN‘^^/2) , ( 2 ) 

where N is an arbitrary parameter that controls the soli- 
ton amplitude and width, ^ = zjr = t/to, = A/Aq, 
^0 = l/(7^o), Zo = 4/\I32\ is the second-order disper¬ 
sion length at ujq, and to is the input pulse duration. 

There are some limitations in using the above NLSE, 
for instance: i) when studying ultrashort pulses, with 
few optical cycles {uoto ^ 1), since the SVEA is no longer 
valid, and ii) investigating polarization effects in birefrin- 
gent waveguides requires to take into account the vector 
nature of the electric field. Alternative, more accurate 
methods of simulating the propagation of electromag¬ 
netic pulses in dielectric media are the finite difference 
time domain (EDTD) method [55ll59] or the recent uni¬ 
directional pulse propagation equation (UPPE) [6Ql l6T] . 
However, accuracy is achieved at the expense of an in¬ 
creased computational effort, leading to a deficiency of 
understanding different physical mechanisms behind the 
pulse dynamics. 

There are also several advantages in using the NLSE: i) 
suitability in generalizing the NLSE to include different 
phenomena such as Raman nonlinearity, self steepening, 
and photoionization effects; and (ii) the possibility of us¬ 
ing well-known analytical techniques such as variational 
perturbation theory to study new nonlinear effects m- 
Based on our experience, the NLSE method can usually 
produce very good qualitative and quantitative results, 
usually accompanied by a deep physical understanding of 
pulse dynamics. In this review, it is our aim to present 
recent and novel SVEA equations that are able to pro¬ 
vide a full understanding of the salient dynamics of pulse 
propagation in gas-filled fibers. 

Intrapulse Raman scattering redshift — In a molecu¬ 
lar medium, a fraction of optical power can be transferred 
from one pulse to another via Raman effect, when the fre¬ 
quency difference between the two pulses matches with 
the vibrational modes of the medium. This effect can oc¬ 
cur within a single pulse, when it is too short and has a 
broad spectrum that exceeds the Raman-frequency shift. 
In this case, the high-frequency (blue) spectral compo¬ 
nents of the pulse continues to amplify the low-frequency 
(red) components during propagation. This amplifica¬ 
tion appears as a redshift of the pulse spectrum, known 
as intrapulse Raman redshift. The NLSE can be modified 
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to study this effect in silica-core fibers m, 

m=2 ^0 

where tr is the Raman coefficient in normalized units. 
To investigate this effect, higher-order dispersion co¬ 
efficients are first neglected {Prn >2 ~ 0) for simplicity. 
For weak Raman nonlinearity, and absence of higher- 
order dispersion coefficients Prn> 2 i the solution of Eq. 
(© can still be assumed to be a fundamental soliton 
that is perturbed by the Raman effect, i.e. 'ip (^, r) = 
N sech [N {t — f (^))] exp [—i5 (^) r], where 5 and r are 
the soliton central frequency and temporal peak that 
change during propagation because of Raman scattering. 
Applying perturbation theory m, the soliton is found 
to be linearly redshifting in the frequency domain with 
rate g = and decelerating in the time domain, 

i.e. S = and f = g(^j^. 


III. RAMAN EFFECT IN GAS-FILLED 
HC-PCFS 

Unlike silica glass, which has a very broad Raman spec¬ 
trum, stimulated Raman scattering processes in gases 
have a narrow Raman-gain spectrum. For this reason, 
gases are characterized by having a very long molecu¬ 
lar coherence relaxation (dephasing) time, of the order 
of hundreds of picoseconds or even more, in compari¬ 
son to the short relaxation time of phonon oscillations 
in silica glass (approx. 32 fs). Within this long re¬ 
laxation stage, the medium can exhibit a highly non- 
instantaneous response to pulsed excitations. Raman 
responses can be manifested in either rotational or vi¬ 
brational modes. Nonlinear interactions between opti¬ 
cal pulses and Raman-active gases have been usually ex¬ 
ploited in the synthesis of subfemtosecond pulses j62f[66] . 

Conventional approaches have been used to enhance 
Raman scattering in gases such as focusing a laser beam 
into a gas confined in a bore-fiber capillary m , or em¬ 
ploying a gas-filled high-finessse Fabry-Perot cavity to 
increase the interaction length [68]. However, Benabid et 
al have exploited the benefits of HC-PCFs and have 
demonstrated Stokes-generation via Raman scattering 
using only few micro joule optical pulses in hydrogen- 
filled HC-PCF [35]. The range of the used energies in 
this experiment were nearly two-order of magnitude less 
than other values reported in prior techniques; demon¬ 
strating the capability of HC-PCFs to enhance nonlinear 
Raman interactions in gases. 

F. Belli et al. have shown an ultrabroadband super con¬ 
tinuum generation spanning from 125 nm in the vacuum- 
UV to 1200 nm in the mid infrared regime, when pump¬ 
ing a H 2 -filled photonic crystal fiber using a 30 fs pulse 
centered at wavelength 805 nm with energy 2.5 /iJ [69] . 
The uniqueness of this work is the extension of supercon¬ 
tinuum generation below 200 nm. Due to the interplay 


between the Kerr effect, the Raman responses of both 
the rotational and vibrational excitations, and the shock 
effect along the fiber with different levels of pulse com¬ 
pression, a dispersive wave at 182 nm on the trailing edge 
of the pulse is emitted that broadens the spectrum into 
the vacuum-UV region. 

Density matrix theory — The dynamics of the Raman 
polarization (also called eoherenee) Pr due to a single 
Raman mode excitation in gases can be determined by 
solving the Bloch equations for an effective two-level sys¬ 
tem [66l [To] . 


dtw + 


re + 1 iai2 


Ti 


{p,2-pl2)E^ 


dt + — - iujR 
J-2 


h 

'i 

P12 = ^ [ai2W + {ail - 0122) P12] E'^, 


(4) 

where aij and pij are the elements of the 2x2 polarizabil¬ 
ity and density matrices, respectively, E (z, t) is the real 
electric field, cjr is the Raman frequency of the transition, 
w = P 22 — pii is the population inversion between the ex¬ 
cited and ground states, P 22 + pii = 1, P 21 = P 125 <^12 = 
0 ^ 21 , Eq is the molecular number density, Ti and T 2 are 
the population and polarization relaxation times, respec¬ 
tively, and h is the reduced Planck’s constant. Solving 
these coupled equations, the Raman polarization is then 
given by Pr = [anpn + a 22 p 22 + ai 2 {pi 2 + pW] EfoE. 
For weak Raman excitation, pn ~ 1 and P 22 ~ 0, i.e. 
the second term in Pr can be neglected, while the first 
term increases the linear refractive index of the medium 
by a fixed amount. 

Using Maxwell and Bloch equations and applying the 
SVEA, one can derive the following set of normalized 
coupled equations that govern pulse propagation in HC- 
PCEs filled by Raman-active gases. 


id^+Yi 


'’"^%ar + |V’P + -Re(pi2) 


‘'O 


m=2 


tVpm\ 




(w 1) to \ I / |2 

OrW H- — -= -4pwlm[pi2) IV^I . 


= 0 , 


dr+^T^-iS 

J-2 


Pl 2 = ipw \tp\ , 


(5) 


where weak Raman excitation is assumed, zr = 
ceo/{ai2 Nqujq) is the nonlinear Raman length, p = 
Po/Pi, Po = ^05 Pi = 2ficeoAeff/((ai2^o), ^ 
and Re and Im represent the real and imaginary parts. 

Eor femtosecond pulses, the relaxation times of the 
population inversion (Ti) and the coherence (T 2 ) can 
be safely neglected, since they are of the order of hun¬ 
dreds of picoseconds or more. Eor instance, Ti 20 
ns and T 2 = 433 ps for excited rotational Raman in 
molecular hydrogen under gas pressure 7 bar at room 
temperature [MlIZI]. We have found also that the pop¬ 
ulation inversion is almost unchanged from its initial 
value for pulses with energies in the order of few pJ, 
i.e. w (r) ~ w(—oo) = —1. The set of the governing 
equations Eq. ® can be reduced to a single generalized 
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FIG. 1. Temporal evolution of an accelerated oscillating Ra¬ 
man polarization with period A = 56.7 fs induced by a prop¬ 
agating fundamental soliton with an amplitude Ni = 1.33, 
a central wavelength 1064 nm, and a FWHM 15 fs in a H 2 - 
filled HC-PCF with a Kagome-lattice cross section, a flat-to- 
flat core diameter 18 /xm, a zero dispersion wavelength 413 
nm, a gas pressure 7 bar, and a rotational Raman frequency 
Ct;R = 17.6 THz. The dashed yellow line represents the tempo¬ 
ral evolution of the soliton that excites the coherence wave. 
The simulation parameters are 7 = 7.07 x 10“® 

^2 — —3425.5 fs^/m, Agff = 134/xm^, ai 2 = 0.8 x 10 “^^ C 
m^/V [69l[73], and to = 11.34 fs. The parameter 7 is calcu¬ 
lated using the nonlinear susceptibility of H 2 m- For these 
parameters, we have found that higher-order dispersion and 
self steepening effects have a weak influence on the soliton 
dynamics. 


nonlinear Schrodinger equation, 

id^i> + i^^drip + + 1^1^^ + -R (r) t/) = 0, (6) 

to ^ 

where R{t) = n sin [5 (r — r')] (t') dr' is the re¬ 
sulting Raman convolution, and k. = is the ratio 

between the Raman and the Kerr nonlinearities. Pump¬ 
ing in the deep anomalous dispersion regime (/32 < 0), 
it is assumed that higher-order dispersion coefficients 
&m >2 can be neglected. For ultrashort pulses with du¬ 
rations to ^ 1/cc^R, sin [(i (r — r')] can be expanded 
around the temporal location of the pulse peak by us¬ 
ing the Taylor expansion. For instance, a fundamental 
soliton with amplitude N and centered at r = 0 will 
induce a Raman contribution in the form of R (r) ~ 
sin ((ir) [1 + tanh (Atr)] at the zeroth-order Taylor 
approximation. This soliton will generate a retarded si¬ 
nusoidal Raman polarization that can impact the dynam¬ 
ics of the other trailing probe pulse lagging behind it. On 
the other hand, for to ^ 1 /^r, ~ 7r (^)I^ with 

7r = i-c the Raman nonlinearity is considered to be 
instantaneous. So, the Raman contribution would induce 
an effective Kerr nonlinearity that is significant, and can 
compete directly with the intrinsic Kerr nonlinearity of 
the gas |69l[72] . 

In the following, we will exploit the long Raman 


coherence induced by an ultrashort pulse in controlling 
pulse dynamics, see Fig. We will study the prop¬ 
agation of two successive pulses separated by a delay 
^ Ti,T 2 in the deep anomalous dispersion regime. The 
two pulses are assumed to have the same frequency, 
hence, they will propagate with the same group velocity, 
and experiencing the same dispersion. The leading 
pulse is an ultrashort strong ‘pump’ pulse ipi with 
to ^ 1/^R. In this case, Eq. ^ can be used to 
determine the pump solution by replacing -0 by 
For weak Raman nonlinearity, the solution of Eq. 
can be assumed to be a fundamental soliton that is 
perturbed by the Raman polarization, i.e. ijji (^, r) = 
Ni sech [Ni (r - ui^ - fi (O] exp [-iQi (C (r - uiC)], 
where Ui = /SuZo/to, Pn is the first-order dispersion 
coefficient of the pump, A^i, and fi are the soliton 
amplitude, central frequency, and temporal location 
of the peak maximum, respectively. Assuming that 
we launch this soliton as a pump with ri (0) = 0, and 
using the variational perturbation method m, we have 
found that this soliton is linearly redshifting in the 
frequency domain with rate gi = ^K7r5‘^csch{7r6/2Ni)^ 
and decelerating in the time domain, i.e. = —gi 
and fi = gi^‘^12. In the case to < I/ciJr, we have 


found that a factor of 


might be used to correct 


the overestimated value of ^i, resulting from using the 
zeroth-order Taylor approximation. The treatments 
of the dynamics of the trailing pulse ‘probe’ ^2 are 
presented in Sec. [Ill A and [III B] when it is a weak long 
pulse and strong ultrashort soliton, respectively [53l [54] . 
The study of the pump-probe dynamics will prove crucial 
in the future understanding of the essential building 
blocks of supercontinuum generation in Raman-active 
gases when excited by ultrashort pulses 


A. Weak probe evolution 


When a second weak probe pulse is sent after the lead¬ 
ing pump soliton, the probe evolution is ruled by the 
following equation: 

+ iu 2 dr ^2 + + 2k.Ni sin {5f) 7/^2 = 0, (7) 


where U 2 = /3i2^oAo, rn = |/32i| / 1/^221, Pij and hj are 
the first and the second order dispersion coefficients of 
the pulse with j = 1,2. Going to the reference frame 
of the leading decelerating soliton, f = r — xxi^ — ^i^^/2, 
and applying a generalized form of the Gagnon- 
Belanger phase transformation m ^2(e,7 = 


(j> {£,, f) exp if (gi^ + ui - U 2 ) + i (gi^ + ui - U 2 f /&gi 
Eq. 0 becomes [53] 


id^cj) = — 


1 

2m 


df(j) + [—2kNi sin {5f) + gif] (j). 


(8) 


This equation is the exact analogue of the time-dependent 
Schrodinger equation of an electron in a periodic crystal 
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FIG. 2. A portion of the absolute eigenstates of a Raman- 
induced temporal periodic crystals with a lattice constant A = 
56.7 fs in the presence of a force with magnitude gi — 0.1408 
in the positive-delay direction. The vertical axis represents 
the corresponding eigenvalues —q. The dotted-dashed line 
is the potential under the applied force. Other simulation 
parameters are similar to that used in Fig. 

in the presence of an external electric field. In Eq. © 
time and space are swapped with respect to the con¬ 
densed matter physics system, as usual in optics, and we 
deal with a spatial-dependent Schrodinger equation of a 
single particle ‘probe’ with mass m in a temporal crys¬ 
tal with a periodic potential U = —2hzNi sin {6f) in the 
presence of a constant force —gi in the positive-delay di¬ 
rection. The leading soliton excites a sinusoidal Raman 
oscillation that forms a periodic structure in the reference 
frame of the soliton, as shown in Fig. Due to soliton 
acceleration induced by the strong spectral redshift, a 
constant force is applied on this structure. Substituting 
4>{i,T) = f (f) exp {iqC), Eq. ^ becomes an eigenvalue 
problem with eigenfunctions /, and eigenvalues —q. The 
modes of this equation are the Wannier functions ITT] 
that can exhibit Bloch oscillations m, intrawell oscilla¬ 
tions eg, and Zener tunneling m due to the applied 
force. 

Wannier-Stark ladder — Consider the propagation of 
an ultrashort soliton with FWHM 15 fs in the deep 
anomalous dispersion regime of a H 2 -filled HC-PCF with 
a Kagome lattice. Exciting the rotational Raman shift 
frequency in the fiber via this soliton will induce a long- 
lived trailing temporal periodic crystal with a lattice con¬ 
stant A = 56.7 fs, see Eig. corresponding to the time 
required by the H 2 molecule to complete one cycle of 
rotation. In the absence of the applied force, the solu¬ 
tions are the Bloch modes, while in the presence of the 
applied force, the periodic potential is tilted, and the 
eigenstates of the system are the Wannier functions por¬ 
trayed as a 2D color plot in Eig. where the horizontal 
axis is the time and the vertical axis is the correspond¬ 
ing eigenvalue. These functions are modified Airy beams 
that have strong or weak oscillating decaying tails. Af¬ 


ter an eigenvalue step giA, the eigenstates are repeated, 
but shifted by A, forming the Wannier-Stark ladder, well- 
known in condensed matter physics. As shown, each po¬ 
tential minimum can allow a single localized state with 
very weak tails. A large number of delocalized modes 
with long and strong tails exist between the localized 
states. 

Bloch oscillations and Zener tunneling — An arbi¬ 
trary weak probe following the soliton will be decom¬ 
posed into the Wannier modes of the periodic temporal 
crystal. Due to beating between similar eigenstates in 
different potential wells, Bloch oscillations arise with a 
period S/gi^ while beating between different eigenstates 
in the same potential minimum can result in intrawell 
oscillations. In our case we did not observe in the sim¬ 
ulations the latter kind of beating, since only a single 
eigenstate is allowed within each well. Interference be¬ 
tween modes lying between different wells are responsible 
for Zener tunneling that allows transitions between differ¬ 
ent wells (or bands). In the absence of the applied force 
{gi = 0)5 the band structure of the periodic medium can 
be constructed by plotting the propagation constants of 
the Bloch modes over the first Brillouin zone [—J/2, J/2], 
as shown in Eig. |^a). Zener tunneling occurs when a 
particle transits from the lowest band to the next-higher 
band. The evolution of a delayed probe in the form of 
the first Bloch mode inside a H 2 -filled HC-PCE under the 
influence of the pump-induced temporal periodic crystal, 
is depicted in Eig. [^b). Portions of the probe are lo¬ 
calized in different temporal wells. Bloch oscillations are 
also shown with period 34.7 cm, which correspond to the 
beating between localized modes in adjacent wells. Af¬ 
ter each half of this period, an accelerated radiation to 
the left due to Zener tunneling is also emitted. Zener 
tunneling is dominant, and Bloch oscillations are weak, 
because the potential wells are relatively far from each 
other. The overlapping between the localized modes are 
small, consistent with the shallowness of the first band 
in the periodic limit (absence of the applied force). 


B. Strong probe evolution 

We now study the case when the delayed pulse is an¬ 
other strong fundamental soliton, rather than a weak 
probe pulse as in [53]. The governing equation for this 
strong ‘probe’ is given by [54| 

id^'lp2 + + |V’2|^V’2 + R2 (r) 1p2 = 0 . ( 9 ) 

Eor weak Raman nonlinearities, the solution of 
this equation is another perturbed fundamental soli- 
ton, i >2 (C, t) = N 2 sech [^2 (t - f 2 (C))] exp [-iQ ,2 (C) r] 
where N 2 ^ and r 2 are the second soliton ampli¬ 
tude, central frequency, and temporal location of the 
peak maximum, respectively. When the soliton dura¬ 
tion ^1/(5, its Raman response function can be approx- 
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(b) Temporal evolution 
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FIG. 3. (a) Bandstructure of the temporal crystal induced 
by the leading ultrashort soliton propagating in the H 2 -filled 
Kagome-lattice HC-PCF with m = 1 in the absence of the 
applied force, (b) Temporal evolution of a weak probe in 
the accelerated periodic temporal crystal. The probe initial 
temporal profile is a Gaussian pulse with FWHM 133.6 fs su¬ 
perimposed on the first Bloch mode of the periodic crystal in 
the absence of the applied force. The contour plot is given in 
a logarithmic scale and truncated at -60 dB. Other simulation 
parameters are similar to that used in Fig. 


imated by using a Taylor expansion as [53] . 

i ?2 (t) « K ^ Ni sin [(5 (r - f;)] {1 + tanh [Ni (r - f;)]} . 
1 = 1,2 

( 10 ) 

The superposition of the induced-Raman effects by the 
two solitons will affect the trailing soliton dynamics. 

Adopting the variational perturbation method to un¬ 
derstand how Raman nonlinearities can affect the pulse 
dynamics m, we have derived a set of coupled govern¬ 
ing equations that determine the evolution of each soliton 
parameters [54], the solutions of which are 


— —91 

Ti = gi C^/2, 


^2 =-52^-52^/ COS [(5 (r2 - Ti)] (11) 

-'''2 Jo 


T2 = 


/ 




where gj = ^K7r6‘^csch{7r6/2Nj). The first (leading) soli¬ 
ton will always linearly redshift in the frequency domain 
with rate gi , and decelerate in the time domain. Whereas 
for the second (trailing) soliton, its dynamics depends on 
two different components: (i) its own (self) component 
that will lead to a linear redshift similar to the leading 
soliton, with rate ^ 2 ; (ii) a cross component represent¬ 
ing the effect of the first soliton on the second soliton. 
The latter component is proportional to the ratio be¬ 
tween their amplitudes and the cosine of the time differ¬ 
ence between them. Since the cosine term varies between 
positive and negative values, the dynamics of the second 


Temporal 

Si^tio^mf^ral 

crystal 

crystal 


\ ' 

Well 

Barrier 


Jc)_ 




FIG. 4. (Golor online). Raman polarization induced by two 
fundamental solitons propagating in a gas-filled HG-PGF in a 
reference frame moving with the leading soliton. White solid 
lines represent the temporal evolution of the two solitons with 
normalized amplitudes Ni = 2.5, N 2 = 1.25 ( corresponding 
to full width half maximum (FWHM) 8 fs, 16 fs, respectively), 
using the analytical prediction Eq. ( [TT| . The second soli¬ 
ton is launched at (a) T 2 (0) = 6.75A,*'^) T 2 (0) = 6.5A, (c) 
T 2 (0) = 5 . 25 A, and (d) T 2 (0) = 5A, where A = 27r/(5, S = 0.5 
(equivalent to a Raman-mode oscillation with period 185 fs, 
such as in deuterium ED), and K, = 0.16. 


soliton can switch back and forth between redshift and 
blueshift in the frequency domain or deceleration and 
acceleration in the time domain. This analytical model 
shows a very good agreement with the numerical model 
provided that the assumption of the soliton durations 
<C 1/d is satisfied. This method can be extended easily 
to the case of more than two solitons. 

Figure]^ shows four special cases of the temporal evo¬ 
lution of the second soliton superimposed on the total 
induced-Raman polarization in a reference frame mov¬ 
ing with the leading soliton f = r — ^i^^/2. A ^2 is chosen 
smaller than Ni so that the cross component is compara¬ 
ble to the self component in Eq. (11). The trailing soliton 


can be treated as a particle in a moving periodic poten¬ 
tial. As we are operating in the anomalous dispersion 
regime, the positive (negative) variation of the refractive 
index represents a potential well (barrier). Thus, the pe¬ 
riodic modulation of the refractive index corresponds to 
a sequence of alternative wells and barriers. Based on 
the initial time delay between the two solitons Ar^, the 
dynamics of the second soliton behaves differently. Also, 
the uniformity of the temporal crystal will be modified 
along the direction of propagation, resulting in a spa- 
tiotemporal modulation of the refractive index, i.e. a 
spatiotemporal crystal Looking at Fig. launching the 
second soliton at (a) the top of a barrier or (b) at the right 
edge of a well, the second soliton will be able to overcome 
the barriers during propagation, and transported across 
the potential by the acting force to the left direction. The 
output spatiotemporal crystals are chirped along the di- 
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rection of propagation in these cases. Interestingly, the 
second soliton in (b) experiences a net maximum self¬ 
frequency blueshift of 9.12 THz = 51.4 nm after 10 cm of 
propagation, before it is redshifted. Launching the sec¬ 
ond soliton at (c) a potential minimum or (d) a left edge 
of a well, the second soliton will not be able to overcome 
the barriers in these cases, so it is trapped inside the well 
and will oscillate indefinitely. The amplitude of oscilla¬ 
tion in (c) is very small, since the initial velocity of the 
soliton in this potential is zero. The soliton will oscillate 
in an asymmetric manner as in (d) due to the modified 
potential beyond the second soliton peak as well as the 
acting force that is opposite to the initial velocity. The 
resulting spatiotemporal crystals have uniform periods 
along the direction of propagation in these cases. 

The temporal evolution of two successive ultrashort 
Gaussian pulses rather than fundamental solitons are de¬ 
picted in Fig. [^for different time delays. The two pulses 
have the same central frequencies and amplitudes, and 
the delay is again within the relaxation coherence time 
of the Raman-active gas. The two pulses will experi¬ 
ence pulse compression and soliton fission processes. The 
‘first’ leading pulse excites the Raman polarization ‘po¬ 
tential’ that will affect the ‘second’ trailing pulse. The 
dynamics of the leading pulse is certainly independent of 
the time delay and will encounter a self-induced Raman 
redshift (deceleration). The trailing pulse dynamics is 
influenced by its self Raman-induced effect as well as the 
cross-accelerated Raman polarization effect of the lead¬ 
ing pulse. In Fig. panel (a) shows the case when the 
self and cross components are working together, result¬ 
ing in a strong delay in comparison to the first pulse. 
The situation when the self and cross components acts 
against each other is figured in panel (b), where initially 
the second pulse is nearly halted since the cross and self 
components cancel each other. When the second pulse 
is launched at one minima of the potential induced by 
the first pulse, the pulse is well-confined during propaga¬ 
tion. Even after the pulse fission, the generated solitons 
are still traveling together, see panel (c). Launching at 
one potential-maxima, the dynamics of a tree-like behav¬ 
ior is obtained as shown in panel (d), where each soliton 
propagates in a different direction. In all these cases, 
the dynamics of each soliton depends on where exactly 
this soliton is born inside the total accelerated periodic 
potentials induced by other preceding solitons. 

The above results will be of fundamental importance 
in the understanding of the building blocks of supercon¬ 
tinuum generation: the multitude of solitons propagating 
in the fiber influence each other in a well-defined way, de¬ 
fined by their intensities and their temporal separations. 
The condensed matter physics analogue effects will take 
place and they will determine the dynamics of each in¬ 
dividual solitons in the supercontinuum process [Ml [75]. 
This theory allows new nonlinear phenomena that are 
impossible to achieve in conventional solid-core optical 
fibers, and opens up new exciting venues for future dis¬ 
coveries. 





t (ps) t (ps) 


FIG. 5. (Color online). Temporal evolution of two successive 
identical Gaussian pulses with profile iVexp(— t^/2 ) in the 
gas-filled HC-PCF described in Fig. [^ with N = 7, and 
FWHM = 25 fs. The first pulse is launched at ti (0) = 0. The 
second pulse is launched at: (a) t 2 (0) = 2A, (b) f 2 (0) = 2.5A, 
(c) t2 (0) = 2.25A, and (d) t2 (0) = 2.75A. 


IV. IONIZATION EFFECT IN GAS-FILLED 
HC-PCFS 


Ionization Models — Photoionization is the physical 
process in which an electron is released and an ion is 
formed due to the interaction of a photon with an atom 
or a molecule. The free electron density Ue is governed 
by the rate equation [ 8 ^ 

ftn 

=W{t){n-ne)-v^e-Prnl, ( 12 ) 

where W is the ionization rate, n is the total density of 
the atoms, and 77 and /3r are the electron attachments, 
and recombination rates, respectively. For pulses with 
pico-second duration or less, both r] and [3r are neg¬ 
ligible. Based on the so-called Keldysh parameter pk 
[ 83] . photoionization can take place by multiphoton ab¬ 
sorption Pk ^ 1 , tunneling ionization px ^ 1 or both 
when pk ^ 1. Several models have been developed to 
determine the dependence of the ionization rate on the 
electric field of the optical pulse, for instance, multi¬ 
photon ionization-based Keldysh-Faisal-Reiss model [84| , 
tunneling based Ammosov- Delone-Krainov method [85] , 
Perelomov, Popov, and Terentev (PPT) hybrid technique 
[ 86 ] . and Yudin-Ivanov model m, which is a modifica¬ 
tion of the PPT technique. It has been shown experi¬ 
mentally that tunneling ionization is dominant over the 
multiphoton ionization in noble gases for pulses with in¬ 
tensities r\j 10 ^^ W/m^ jHHUM], which is our case of study. 
A review of these processes and models is found in Ref. 
m- In the tunneling regime, the time-averaged ioniza¬ 
tion rate is given by [82] 

W{I) = d{lHliy/^ exp[-b{lH/I)^^% (13) 










where d = 4So [ 3 / 7 r]i /2 [Ui/UhV/\ b = 2/Z [Ui/Uh?/\ 
(5o = 4.1 X 10^^ Hz is the characteristic atomic frequency, 
[//is the ionization energy of the gas, Uh ~ 13.6 eV is 
the ionization energy of atomic hydrogen, Ih = 3.6x10^^ 
W/cm^ and / = is the laser pulse intensity. This 
model is adequate for the analysis that is based on the 
evolution of the complex envelope of the pulse. Un¬ 
fortunately, all these models are not straightforwardly 
amenable to analytical manipulation, because of the com¬ 
plex dependence on the pulse intensity. Equation (13) 


predicts an ionization rate that is exponential-like for 
pulse intensities above a threshold value Uh [ISl [9T] . Any 
pulse with an intensity / Uh will suffer a strong ioniza¬ 
tion loss due to the absorption of photons in the plasma 
generation process. This limits the operating regime to 
near the ionization threshold Uh where the ionization loss 
is drastically reduced. A model of the ionization rate with 
linear dependence on the pulse intensity can be developed 
using the first-order Taylor series of Eq. (13) [49] 


W^d(/-Uh)0(/-Uh), 


(14) 


where d is a constant that is chosen to reproduce the 
physically observed value of the ionization threshold, and 
0 is a Heaviside function, introduced to cut the ionization 
rate to zero below the value of Uh- 

Photoionization results in decreasing the refractive in¬ 
dex of the medium by a factor proportional to the square 
of the plasma frequency cjp, also it attenuates the pulse 
amplitude due to photon absorption. To include these 
effects Eq. 0 can be modified as [49] 


id,A+Y\ —,PmdrA+-i\A\^A---^+i 


.AefiUi 


m=2 


m\ 


2cJoc 2\A\ 


dtUe = 0, 

(15) 

where c is the speed of light, and Aeff is the effective 
area. Using the split-step Eourier method, the pulse am¬ 
plitude can be determined at each propagation step after 
computing the free electron density Uq. After neglect¬ 
ing the electron attachments and recombination effects, 
Eq. (12) can be directly integrated analytically, and can 
then be substituted in Eq. ( 01 ) to have a single gen¬ 
eralized NLSE to describe pulse propagation in an ion¬ 
izing medium. Introducing the following rescalings and 
redefinitions: = e^n/ [eorue] is the maximum plasma 

frequency, 0 = ^koZo [cjp/cjo]^, [^t/^o]^, 

O' = d-U/[Aefr7K^Q], and k = Ui a eo rUe uj^/[ko e‘^] 


In 


this case, Eqs. ( T3p5 ) become [9T] 


id^2p - 




m=2 ^ 


+ + = 0, (16) 


where 0 = 


1 — exp 




(a = (0T - 0) A|Ur0A|Un and A|Ur = “ IV'lth’ 


A. Short pulse evolution 

In order to extract useful analytical information from 
Eq. (16), further simplifications are necessary. Eirst, 


we assume operating in the deep anomalous regime. Eor 
pulses with maximum intensities just above the ioniza¬ 
tion threshold, also called floating pulses m, the ion¬ 
ization loss is not large and can be neglected as a first 
approximation. Eor such pulses, only a small portion of 
energy above the threshold intensity contributes to the 
creation of free electrons. Eurthermore, the effect of the 
0-function can be approximately determined via multi¬ 
plying the cross-section parameter cr by a factor 5 that 
represents the ratio between the pulse energy contribut¬ 
ing to plasma formation (the portion above the ioniza¬ 
tion threshold) and the total energy of the pulse. In 
this case, ionization can be treated as a perturbation 
of the solution of the NLSE, which is the fundamen¬ 
tal soliton. The solution of Eq. (16) can be written 
as U (^5 = ^(0 sech where r 

is the temporal location of the soliton peak, and 6 is 
the pulse central-frequency shift. Applying the pertur¬ 
bation theory [T6|, 6{^) = —and f (^) = g^ where 
g = — (2/3)5cr0TA'^ [49]. This shows that photoioniza¬ 
tion should lead to an absolutely remarkable soliton self¬ 
frequency hlueshift. This blueshift is accompanied by a 
constant acceleration of the pulse in the time domain - 
opposite to the Raman effect, which produces pulse de¬ 
celeration [T5] . 

Including the effects of both the photoionization loss 
and the Heaviside function, the perturbation theory re¬ 
sults in a set of two coupled differential equation that 
governs the spatial evolution of the soliton amplitude and 
central frequency. 


= —2k,(I)t (Atanh'd — |U| 

^ = cr0T ^ tanh^'d + 

0 (^ o 


T) 


('dsech^'d — tanh'd) 


(17) 

where 'd = AT, and T ^ 1/Asech [|U|th/A] is the 
temporal position where the pulse amplitude exceeds the 
ionization threshold at ^ = 0 [91]. Solving these equa¬ 
tions numerically as shown in Eig. we found that 
pulses with initially large intensities (Aq > |Ulth) 
experience a boosted self-frequency blueshift. However, 
the ionization loss suppresses the soliton intensity after a 
short propagation distance to the floating-soliton regime, 
where the soliton can propagate for a long propagation 
distance with a limited blueshift and negligible loss. The 
maximum frequency shift is achieved when the soliton 
intensity falls below the photoionization threshold. 

To study the full dynamics of pulse propagation in 
fibers filled by an ionizing gas, Eq. (16) should be solved 


Ith 


= Uh^eff- 


numerically via the split-step Eourier method. The tem¬ 
poral and spectral evolution of a pulse, with an initial 
temporal profile Asech (r) and intensity less than the 
ionization threshold, are depicted in the panels (a,b) of 
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FIG. 6. (Color online). The spatial dependence of the soli- 
ton amplitude and frequency shift during a photoionization 
process for different initial pulse amplitudes. 


Fig. respectively. Panel (c) shows the variation of the 
ionization fraction along the fiber. The pulse is pumped 
in the deep anomalous-dispersion regime of the fiber, and 
it undergoes self-compression. When the pulse intensity 
exceeds the threshold value, a certain amount of plasma 
is generated due to gas ionization, and a fundamental 
soliton is ejected from the input pulse. The soliton cen¬ 
tral frequency continues to shift towards the blue-side due 
to the energy received from the generated plasma. How¬ 
ever, because of the concurrent ionization loss, the soli¬ 
ton intensity gradually attenuated to the regime where 
Such pulses, the floating solitons, can prop¬ 
agate for considerably long distances with minimal atten¬ 
uation and limited blueshift. When the soliton intensity 
goes below the ionization threshold, the blueshift pro¬ 
cess is ceased. A second ionization event accompanied 
by a second-soliton emission can take place by further 
self-compression of the input pulse based on its initial 
intensity. At the end, a train of floating solitons are gen¬ 
erated. A clear representation for the pulse dynamics in 
the presence of plasma is shown in Fig. where the 
temporal profile of the pulse intensity is plotted at 
selected positions inside the fiber. 

Long-range non-local soliton forces and clustering — 

A non-local interaction between two successive soli¬ 
tons have been found when their temporal separation 
is shorter than the recombination time [49], similar to 
the nonlinear interactions in Raman-active gases within 
the molecular coherence relaxation time presented in Sec. 
III. The leading soliton and its induced non-vanishing 
electron-density tail both co-propagate and accelerate to¬ 
wards the negative-delay direction. Within the recom¬ 
bination time, the trailing soliton will be affected by a 
force in the opposite direction due to the accelerated 
long electron-density tail. The acceleration of the trailing 
soliton acquires an exponentially decaying dependence on 
the amplitude of the leading soliton. These dynamics are 
featured in Fig. [^ that shows the temporal and spectral 


Temporal evolution 







FIG. 7. (Golor online). Temporal (a) and spectral (b) evolu¬ 
tion of an ultrashort pulse with in an Ar-filled HG-PGF. The 
temporal profile of the input pulse is N sechr, with iV = 8, 
to = 50 fs. The gas pressure is 5 bar. Gontour plots are 
given in a logarithmic scale, (c) Spatial dependence of the 
ionization fraction along the fiber. 


dependence on the soliton parameter N assuming that 
the input pulse is N sechr. The scenario is as follows: As 
long as the intensity of the first-emitted soliton is above 
the threshold intensity, it prevents the ejection of a sec¬ 
ond soliton due to the presence of the opponent force 
of the non-local interaction. When the ionization loss 
ends the blueshifting process of the leading soliton, the 
trailing soliton can be ejected and recovers its expected 
acceleration and blueshift. This allows the second soli¬ 
ton to catch up and cluster with the first soliton. In 
addition, the spectrum of the two solitons start to over¬ 
lap and form spectral clustering. Similarly, when the 
first two solitons are very close to each other, their in¬ 
duced electron-density tail applies a combined force on 
the third soliton. The evolution of the cross-frequency- 
resolved optical gating (XFROG) spectrograms for pulses 
with different initial intensities are depicted in the panels 
of Fig. 10, where (a) represents the input pulse; (b) and 
(c) shows the emission of the first and second solitons, 
respectively; and (d) depicts the temporal and spectral 
clustering of the first two solitons and the emission of a 
third soliton. 
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(a) Initial pulse 


^ = 0 


500 


(b) First self-compression 


(c) First-soliton ejection 


(d) Second self-compression 


C = 0.14 


(e) Second-soliton ejection 



(f) Floating solitons 




FIG. 8. (Color online). Intensity profile of an optical pulse 
in the time domain at different positions, (f, inside an Ar- 
filled HC-PCF. The dashed red line represents the threshold 
intensity. The simulation parameters are similar to Fig. 
Each panel is titled by its main feature. Insets are enclosed 
in panels for better view and more details. 


FIG. 10. (Color online). XFROG spectrograms for pulses 
with selected soliton parameter N in an increasing order. The 
simulation parameters are similar to Fig. (a) = 3, ^ = 0. 

(b) iV = 5, ^ = 1/4. (c) iV = 8, ^ = 1/4. (d) N = 9.2, ^ = 
1/4. Each panel is titled by its main feature. White arrows 
show the movement of the solitons. 


Temporal evolution 



Spectral evolution 



Soliton order A/ 


EIG. 9. (Color online). Temporal (a) and spectral (b) outputs 
of a pulse with temporal profile N sech r after propagating 
inside an Ar-filled HC-PCE with length ^ = 1/4 versus the 
soliton parameter N. 


B. Long-pulse evolution 

It has been shown recently that when the gas is excited 
by relatively long pulses with ionizing intensities, new 
kinds of self-phase modulation (SPM) and modulational 
instability (MI) can emerge during propagation. More¬ 
over, after the initial stage of instability is over, a ‘shower’ 
of hundreds of solitons, each undergoing an ionization- 
induced self-frequency blueshift, pushes the supercontin¬ 
uum spectrum towards shorter and shorter wavelengths. 
Such a blueshifting plasma-induced continuum has some 


similarities with the redshifting Raman-induced contin¬ 
uum driven by the Raman self-frequency shift in con¬ 
ventional solid-core fibers [92H94] . although the physical 
processes involved are dramatically different. 


Asymmetrical SPM — Ionization-induced SPM can 
be studied analytically using Eq. (16), in the case of 
small dispersion and long input pulse durations the non¬ 
linearity initially dominates over the group-velocity dis¬ 
persion (GVD). Also by temporarily neglecting the losses 
(which do not change the qualitative picture, but only 
saturate the SPM spectrum after a certain distance), and 
assuming Gaussian pulse excitation, close forms of the 
spatial dependence of the mean frequency (ft) and the 


standard deviation (Afl)^ of the output spectrum can 
be derived m- For different values of 77 = which 


depict the spatial dependence of ( 11 ) and (AH)^ along 
the fiber. For 77 = 0, which corresponds to the absence 
of ionization, the mean frequency is always zero during 
propagation due to the well-known symmetric spectral 
broadening due to conventional SPM [16]. As 77 increases, 
the plasma starts to build up inside the fiber. The mean 
frequency moves linearly towards the blue-side of the 
spectrum along the fiber due to the ionization-induced 
phase-modulation. This induces a strong and extremely 
asymmetric SPM, imbalanced towards the blue part of 
the spectrum. In a real case, the spectral broadening 
process is certainly limited by the unavoidable ionization 
and fiber losses. 


measures the ionization strength, panels (a,b) in Fig. 
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Plasma-induced modulational instability — After the 
initial SPM stage is over, the interplay between nonlin¬ 
ear and dispersive effects can lead to an instability that 
modulates the temporal profile of the pulse, creating new 
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FIG. 11. (Color online). Spatial dependence of (a) the mean 
frequency (Q) and (b) the frequency standard deviation AQ of 
a Gaussian pulse exp (—t^/2to) with tq = 2 . The temporal 
position T at which the pulse intensity can initiate photoion¬ 
ization is assumed to be equal to tq. AQq is the spectral 
width at ^ = 0. 


spectral sidebands referred to as modulational instabil¬ 
ity (MI) [95H97| . Starting from Eq. (16), MI due to 
the photoionization nonlinearity can be investigated by 
using the standard approach described in [T6[ |95] . The 
Kerr-induced MI occurs only in the anomalous dispersion 
regime over a defined bandwidth |16]. However, the pres¬ 
ence of the photoionization process induces an unusual 
instability that can exist in both normal and anomalous 
dispersion regimes, and for any frequeney. The spectral 
dependence of the gain of these instabilities on different 
peak powers is shown in Fig. 12 for (a) anomalous and 
(b) normal dispersion regimes, where the physical pow¬ 
ers are normalized to the threshold ionization power, i.e., 
IV^Ith ~ When the normalized input power 
we have the traditional side-lobes, which exist uniquely 
in the anomalous dispersion regime, due to the Kerr- 
nonlinearity. However when ipQ > photoionization- 

induced instability generates unbounded side-lobes with 
slowly-decaying tails. A similar situation occurs in the 
normal regime, however, the gain is slightly lower due to 
the absence of the Kerr contribution. For this case there 
are no instabilities below the threshold power since no 
plasma is generated. 

Propagation of a long Gaussian pulse inside an anoma¬ 
lous dispersive gas-filled HC PCF is portrayed in the pan¬ 
els (a,b) of Fig. 13, obtained by simulating Eq. (16). 
The first stage of propagation shows asymmetric spectral 
broadening towards the blue due to ionization-induced 
SPM. Immediately after the SPM stage, dispersion starts 
to play a role, and due to the combined Kerr and ion¬ 
ization Mis, broad and slowly decaying side lobes are 
generated and amplified quickly. In the third and fi¬ 
nal stage in the propagation, strongly blueshifted soli- 
tons are emitted. Further insight into the dynamics can 


FIG. 12. (Golor online). MI spectral gain profile versus nor¬ 
malized input peak power in (a) anomalous and (b) normal 
dispersion regimes with equal magnitude of GVD. For the 
anomalous regime, a GW with A = 1064 nm is launched into 
a Kagome Ar-filled HG-PGF with core diameter 20 /xm, gas 
pressure of 1 bar, and (32 — —2.8 ps^/km. In this hber, thresh¬ 
old ionization power is ~ 105 MW at room temperature. For 
the normal regime simulation, we assume hypothetically a 
gas-hlled fiber with a ^2 of the same magnitude as in (a), but 
with a positive sign. 


be obtained from Fig. 14, where the evolution of the 
XFROG spectrograms of the pulse at different positions 
along the fiber is shown. The pulse is initially asymmet¬ 
rically chirped in the center of the pulse towards high fre¬ 
quencies, [Fig. (Hb)], due to the higher plasma density 
created at the peak intensities. At the same time two im¬ 
balanced ionization-induced MI sidebands appear in the 
pulse spectrum [Fig. [l^c)]. MI facilitates the formation 
of many solitons. In less than half a meter of propa¬ 
gation the initial pulse disintegrates into a ‘shower’ of 
solitary waves, see Fig. Eld), each undergoing a strong 
self-frequency blueshift induced by the intrapulse pho¬ 
toionization. 


V. CONCLUSIONS AND FINAL REMARKS 

In this review, we have presented very recent re¬ 
sults concerning soliton dynamics in gas-filled HC- 
PCFs. We have divided gases into two main cate¬ 
gories: Raman-active (molecular gases) and Raman- 
inactive (monoatomic gases). The former kind of gases, 
such as molecular hydrogen, are characterized by long 
molecular coherence in comparison to silica glass. This 
results in highly non-instantaneous interactions that can 
be detected by launching a probe pulse, delayed from 
the main pump pulse. For a weak probe, the problem 
is reduced to the motion of a quantum particle in a 
periodic ‘temporal’ crystal subject to an external force. 
Phenomena related to condensed matter physics such as 
Wannier-stark ladder, Bloch oscillations, and Zener tun- 
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neling are predicted to occur. However, if the probe is an¬ 
other ultrashort intense soliton, phenomena such as soli- 
ton oscillations and transport have been shown to occur. 
Moreover, in this case the temporal crystal is upgraded 
to a spatiotemporal one with a uniform or chirped spa¬ 
tial period. In Raman-inactive gases, such as argon, we 
have investigated the effect of photoionization effects on 
the evolution of short and long pulses. Unique phenom¬ 
ena such as soliton self-frequency blueshift, asymmetrical 
self-phase modulation, universal modulation instability, 
and shower of blueshifted solitons have been studied and 
elaborated. These fresh theoretical results, supported 
by recent experiments, can pave the way for the manip¬ 
ulation and control of the pulse dynamics in PCFs for 
demonstrating completely novel optical devices. 


FIG. 13. (a) Spectral and (b) temporal evolution of a long 
Gaussian pulse propagating in an Ar-filled HG-PGF with a 
gas pressure 1 bar, and a hexagonal core-diameter 20 /xm. 
The input pulse has a peak power 200 MW and a duration 
0.707 ps 1/e-intensity half-width. Gontour plots are given in 
a logarithmic scale and truncated at -50 dB. 
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FIG. 14. (Golor online). XFROG spectrograms for the propa¬ 
gation through an Ar-filled HG-PGF. The simulation param¬ 
eters are the same as in Fig. The reference pulse is a 

Gaussian with FWHM 50 fs. (a) Input pulse, (b) SPM fre¬ 
quency chirping, (c) Modulational instability, (d) Pulse dis¬ 
integration into multiple blueshifting solitons. 
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